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ABSTRACT 

The bunching of giant planets at a distance of several stellar radii may be explained by the disruption 
of the inner part of the disk by the magnetosphere of the star during the T Tauri stage of evolution. 
The rotating magnetic field of the star gives rise to a low density magnetospheric gap where stellar 
migration is strongly suppressed. We performed full 3D magnetohydrodynamic simulations of the disk- 
magnetosphere interaction and examined conditions for which the magnetospheric gap is "empty" , by 
changing the misalignment angle between magnetic and rotational axes of the star, 0, and by lowering 
the adiabatic index 7, which mocks up the effect of heat conductivity and cooling. Our simulations show 
that for a wide range of plausible conditions the gap is essentially empty. However, in the case of large 
misalignment angles O, part of the funnel stream is located in the equatorial plane and the gap is not 
empty. Furthermore, if the adiabatic index is small (7 ~ 1.1) and the rotational and magnetic axes 
are almost aligned, then matter penetrates through the magnetosphere due to 3D instabilities forming 
high-density equatorial funnels. For these two limits there is appreciable matter density in the equatorial 
plane of the disk so that a planet may migrate into the star. 



1. INTRODUCTION 

More than 170 giant planets have been discov- 
ered around solar-type stars (see catalog of plan- 
ets in http://vo.obspm.fr/exoplanetes/encyclo/catalog- 
main.php). About 24% of them are located very close 
to the star, at r ^ 0.1AU (see also reviews by Marcy 
et al. 2003; Papaloizou & Terquem 2006). There is a 
prominent peak in the radial distribution of the planets at 
r ~ 0.04 - 0.05 AU ps (7 - 10)i? Q which corresponds to 
periods of 3 days (see Figure 1). 

The radial distribution and other properties of the giant 
planets have stimulated work on models of planet forma- 
tion and migration. According to the presently favored 
interpretation, planets form far away from the star either 
through core accretion (Mizuno 1980; Pollack et al. 1996), 
or through instabilities in the disk (Boss 2001). Subse- 
quently they migrate inward due to their gravitational in- 
teraction with the disk (Lin & Paploizou 1986; Lin, Bo- 
denheimer & Richardson 1996; Ward 1997; Nelson & Pa- 
paloizou 2003). 

Some planets may migrate close to the star where the 
disk properties are strongly influenced by the star and/or 
the star's magnetic field. In particular, the inner regions 
of the disk may be truncated as a result of heating by the 
star (e.g., Kuchner & Lecar 2002) and photoevaporation 
of the region of the disk <; 1 AU (Matsuyama, Johnstone 
& Murray 2003) . Or, the inner part of the disk may be dis- 
rupted by the strong magnetic field of the protostar (Lin 
et al. 1996). 

There are strong observational arguments that young 
solar-type stars (T Tauri type stars) have strong magnetic 
fields. They are thousands of times brighter than the Sun 



in the X-ray (see review by Feigelson & Montmerle 1999) 
which is a sign of their high magnetic activity. In a number 
of cases, direct measurements of the Zeeman broadening 
imply a B field (averaged over the surface of the star) of 
the order of several kilo-Gauss, which is much stronger 
than that the average field on the surface of the Sun (e.g., 
Basri et al. 1992; Johns-Krull et al. 1999). The strongest 
magnetic fields are probably associated with the multipo- 
lar component (e.g. Safier 1998; Johns-Krull ct al. 1999; 
Smirnov et al. 2003). However, a significant dipole com- 
ponent is also expected. It gives many observational signs 
of magnetospheric accretion in the T Tauri stage (see re- 
view by Bouvier et al. 2006) and also is observed in some 
circular polarization measurements (e.g., Valenti & Johns- 
Krull 2004; Symington et al. 2005). 

A sufficiently strong dipole magnetic field will truncate 
the disk at a distance of several stellar radii forming a low- 
density magnetospheric gap as shown in Figure 2. The mi- 
gration rate of planets in this gap will be greatly reduced. 

This work analyzes the conditions where a low-density 
gap exists between the inner edge of the disk and the star's 
surface. Our analysis is based on 3D MHD simulations. In 
particular, we investigate properties of the gaps for differ- 
ent misalignment angles O between the rotation axis of the 
star i?„ and its magnetic moment fi. We investigate con- 
ditions where part of the funnel stream is in the equatorial 
plane close to the star. This may give an appreciable mat- 
ter density in the magnetospheric gap. Also we consider 
the possibility of direct equatorial accretion from the disk 
to the star due to 3D instabilities. 



2. WHERE DOES THE PLANET MIGRATION STOP? 
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Fig. 1. — Distribution of extrasolar planets in the vicinity of the star. 




Fig. 2. — Sketch of an accretion disk which is disrupted by the star's dipole magnetic field. The rate of migration of a 
planet is greatly slowed once it enters the gap. 




Fig. 3. — Radial distribution of the density (solid line) and the angular velocity of the disk lu^ — v^/r (dashed line) in 
the vicinity of the star for a misalignment angle of = 30°. Thin dashed line shows Keplerian angular velocity lok- 



The accretion disk is disrupted at the Alfven radius 
where the dynamic pressure of the disk matter is com- 
parable to the magnetic pressure in the star's dipole field, 



r A = \p? I (MVGM)] 2 / 7 
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where M . 8 = M*/O.8M , R 2 . 5 = R*/2.5 R e , and B 3 = 
i?*/10 3 G are the normalized mass, radius, and magnetic 



Fig. 4. — Result of 3D simulations of disk accretion to a rotating star with a dipole moment misaligned with the star's 
rotation axis f2„ by 9 = 30°. The color background shows the density distribution in the equatorial region. Density 
varies from p as 0.003 (blue color) to p m 1 (red color). The red lines are magnetic field lines. The black arrow shows the 
direction of the magnetic moment fi. 



field of the protostar, and M_ 7 = Af/10" 7 M Q /yr is the 
accretion rate in the disk (Ghosh & Lamb 1979; Camen- 
zind 1990; Konigl 1991). For the typical parameters of 
T Tauri stars used in this formula, this radius coincides 
approximately with the peak of the distribution shown at 
Figure 1. Thus this peak may result from the greatly re- 
duced rate of migration inside the magnetospheric gap. 

Numerical simulations show that the disk is disrupted 
at the distance r « ta, where the plasma is lifted out of 
the disk plane by the vertical pressure force and it then 
flows along the star's dipole field lines in a funnel flow 
(Romanova et al. 2002; 2003; 2004). As a consequence, 
the density of matter in the equatorial plane is greatly re- 
duced for r < ta- Figure 3 shows the equatorial density 
distribution obtained from our 3D simulations. The den- 
sity is large in the disk, and it often increases as ta is 
approached. However, for r < ta the density drops by a 
factor of ~ 100 — 300 in the magnetically-dominated mag- 
netosphere. A protoplanet which migrates inward to radii 
< ta enters a region of greatly reduced density. 

For typical conditions planets migrate inward as a re- 
sult of interaction of the planet with the disk matter. The 
planet loses part of its orbital angular momentum by over- 
taking collisions with the disk outside its orbit and it gains 



a smaller part by overtaking collisions of the disk matter 
inside its orbit. The rate of migration, or radial speed, 
V pr , depends on a number of parameters such as mass of 
the planet, M p , surface density of the disk, S, viscosity 
in the disk, v, and others, and also it is different in cases 
when a planet opens a gap in the disk or not. 

If the planet's mass is relatively small (M p <; 10M®) it 
does not open a gap in the accretion disk. The migration in 
this case is referred to as "type I," and the planet's inward 
drift speed is V pr cx — M p (Sr 2 ) (Ward 1997; Papaloizou & 
Terquem 2006). 

Planets of sufficiently large mass open a gap in the disk 
of width of the order of the disk thickness. The migration 
in this case is referred to as "type II," and it tends to "lock" 
the planet's migration to that of the disk matter if the local 
disk mass, Md — 47rr 2 E , is larger than the planets mass 
M p . The disk matter moves inward with a radial speed 
v r = —3v/(2r), where v = ccc^/uk is the usual Shakura- 
Sunyaev turbulent viscosity with a — 10~ 3 — 10~ 2 . How- 
ever, if the local disk mass is small compared with the 
planet's mass, then the planets migration is slower than 
that of the disk matter. The angular momentum lost by 
the planet dJ p /dt = M p vkV pt /2 in a second is equal to 
the angular momentum transported outward by the vis- 




Fig. 5. — The top panels show the density distributions and sample magnetic field lines in the f2* — \x plane for different 
misalignment angles 9 = 30°, 45°, 60°. The bottom panels show the equatorial distribution of density. 




Fig. 6. — Result of simulations for 7 = 1.1 for 9 = 5° after P — 12 rotations. The left panel shows the density distribution 
and sample field lines in the — fi plane. The middle panel shows the density distribution in the equatorial plane, and 
the right panel shows a 3D view of the density levels and sample field lines. 



cous stress in the disk, dJd/dt = Mtvk, where Vk = 
{GM /r) 1 / 2 is the Keplerian velocity and M mass accretion 
rate of the disk. This gives a migration speed of the planet 
V pr = — (Md/Mp)\v r \. We can write M = 27rr£|u r | so that 
M d = 47rr 2 E « 1.3 x 10 29 g(r/0.1AU) 3 / 2 (M/10" 7 M Q /yr) 
for h/r — 0.1 and a — 1CP 3 . 

For Md < M p ,the time scale for the planet's migration 
at r is simply 



M d iv 



which is independent of r and inversely proportional 



to both M and a. For a Jupiter mass planet, r p w 
4700yr(l(T 7 Af Q /yr/M) for a = 1(T 3 and h/r = 0.1. 
However, equation(2) does not include the influence of the 
star's magnetic field. Inside the magnetospheric gap the 
matter density is reduced by a factor J> 10 2 . The migra- 
tion time will be increased by a corresponding factor. 

Figure 3 also shows that the angular velocity of the disk 
plasma in the equatorial plane Ud within the magneto- 
spheric gap is much lower than the Keplerian angular ve- 
locity of the planet lo k for stars with periods P > 2 days. 
Thus, the disk matter passing close by the planet will have 
a large relative velocity t{ojk ~~ ^d)- The relative veloc- 
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ity is larger than in the non-magnetic case by a factor 
F ~ (r/h)(u>K — ^dj/^K ^ 1 f° r the type II migration. 
One can readily show that the rate of exchange of angular 
momentum between the planet and the disk is reduced by 
a factor 1/F 2 <C 1. Thus the migration time inside the 
magnctosphcric gap will be increased further by a factor 
F 2 . We estimate F 2 ;> 10. Once the planet is inside the 
magnctospheric gap it may continue to migrate slowly in- 
ward owing to resonances it has with the disk matter at 
larger radii but a treatment of this is beyond the scope of 
the present work. 

Thus, in both cases the migration will be greatly re- 
duced if the density in magnetospheric gap is much lower 
than in the disk. Conditions may be changed however if 
the magnetic axis \x is misaligned relative to the rotational 
axis of the star fl* by an angle 6. For high misalignment 
angles 6, the funnel streams may be partially located in 
the equatorial region (Romanova et al. 2003). From the 
other side, even in the aligned case, some matter may ac- 
crete to the star in the equatorial plane due to instabilities 
(e.g., Arons & Lea 1976). We performed 3D simulations 
to investigate both factors. 

3. 3D SIMULATIONS OF THE DISK-MAGNETOSPHERE 
INTERACTION AND MAGNETOSPHERIC GAPS 

We performed a set of 3D simulations using our code 
based on the "cubed sphere" grid (Koldoba et al. 2002) 
with the main goal of analyzing the density distribution in 
the magnetospheric gap. Simulations were set up in a way 
similar to those of Romanova et al. (2003, 2004). Namely, 
quasi-stationary initial conditions were used which per- 
mitted slow viscous accretion from the disk to a star. An 
a— viscosity was incorporated to the code with typical val- 
ues of a— parameter: a = 0.02 and 0.04. The magnetic 
axis fx is misaligned relative to rotational axis of the star 
fl* by an angle 6. The rotational axis of the star coincides 
with that of the disk. 

Simulations were done for parameters typical for T Tauri 
type stars: M* = 0.8M Q , R* = 2.5R , = 10 3 G, 
M w 3 x 10 _8 M Q /yr. Compared to our previous 3D runs, 
we changed parameters so as to increase the size of mag- 
netospheric gap to ta — (4 — 5)i?* versus ta = (2 — 3)ii* 
in our previous papers. Figure 4 shows the magneto- 
spheric gap in a test case with an even larger magneto- 
sphere ta = (6 — 7)ii*. One can see that the low-density 
magnctospheric gap can be quite large. One can also see 
that inside the magnctospheric radius (which corresponds 
approximately to the edge of the disk), the magnetic field 
lines are closed, while outside of this radius they are carried 
by the matter of the disk or corona. Our simulations show 
that the density inside the magnetospheric gap is about 
100 — 300 times smaller than the density in the nearby 
disk. 

Magnetospheric Gap at Different O. We performed 
3D simulations for different misalignment angles O from 
O = 0° to 90°. We investigated the magnctospheric gaps 
in the equatorial plane. Simulations have shown that mat- 
ter flow is different at small and large misalignment an- 
gles. For angles 8 <; 45°, matter flows to the star along 
funnel streams which are above and below the equatorial 
plane. Thus within the magnetospheric gap r < ta the 
matter density in the equatorial plane is greatly reduced 



(see Figure 5, left two panels). For larger angles, matter 
also accretes to a star through the funnel streams. How- 
ever, part of the funnel streams is located in the equatorial 
plane and the magnetospheric gap is not empty (Figure 5, 
right panel). Thus at large O, the planets orbiting in the 
equatorial plane will interact with the dense gas of the 
streams and may continue to migrate inward to the star. 
Accretion through Equatorial Funnels at Low 7. 
There is another possible reason why the magnetospheric 
gap may have some matter density. There are possible in- 
stabilities which may lead to the direct accretion of matter 
through the magnetosphcrc in the equatorial plane. To 
investigate such instabilities, we took the almost aligned 
case, = 5°, and decreased the adiabatic index from 
7 = 5/3 to 7 = 1.1. The adiabatic index may be sig- 
nificantly lower than its ideal value in the case of high 
electron heat conductivity which may occur in a highly 
ionized plasma. In our simulations, the low value of 7 acts 
to give a low temperature in the disk and the funnel flow. 
We observed that matter partially accretes in the equato- 
rial plane. Figure 6 shows that matter accreted through 
funnels which are located inside the magnetosphere. They 
penetrate inwards through the Rayleigh- Taylor type in- 
stability (e.g. Arons & Lea 1976) up to some distance 
?T, and then form regular funnel streams along the field 
lines. The distance of penetration depends on the ratio 
ta/R*- At relatively small values ta/R*, the equatorial 
funnels may penetrate almost to the surface of the star, 
as shown in Figure 6. At larger values of ta/R*, the fun- 
nels move inward only part of the way. Thus, in the case 
of a weak magnetic field and/or high accretion rate, the 
magnetospheric radius ta ~ (1 — 3)i?*, and a small adia- 
batic indexes 7, the gap will not be empty and planets will 
continue to migrate inward unless the tidal interaction or 
some other force will prevent them against falling to the 
star. In the opposite case of a larger magnetosphere, the 
equatorial funnels will occupy only a part of the gap and 
planet may survive longer inside the innermost gap. 

4. CONCLUSIONS 

We conclude that for typical parameters of solar-type 
protostars a very low-density magnetospheric gap forms. 
We show that the rate of inward migration of protoplan- 
ets is greatly reduced within this gap. This gap may ex- 
plain the observed accumulation of planets at the distances 
<~ (0.04 — 0.05) AU from the star. Note however that there 
may be significant matter density within ta in cases of high 
misalignment angles O ;> 45° where the part of the funnel 
stream is located in the equatorial plane. Further, for low 
values of the adiabatic index, 7 = 1.1, and small misalign- 
ment angles, O <; 5°, matter may accrete in the equatorial 
plane due to instabilities. However, the azimuthal veloc- 
ity of the matter within ta differs substantially from the 
Keplerian angular velocity of the planet. This difference 
in the velocities acts to slow the planet's migration. 

We thank Drs. G.V. Ustyugova and A.V. Koldoba for 
valuable contributions to our MHD simulation codes. Also 
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